Second sound dipole mode in a partially Bose-Einstein condensed gas 
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We study the second sound dipole mode in a partially Bose-Einstein condensed gas. This mode 
is excited by spatially separating and releasing the center-of-mass of the Bose-Einstein condensate 
(BEC) with respect to the thermal cloud, after which the equilibration is observed. The oscillation 
frequency and the damping rate of this mode is studied for different harmonic confinements and 
temperatures. The measured damping rates close to the collisionless regime are found to be in 
good agreement with Landau damping. For increasing hydrodynamicity of the cloud we observe an 
increase of the damping. 



PACS numbers: 03.75.Kk, 47.37.+q, 67.85. De 

In 1938 Rapitza, and independently Allen and Mis- 
ener, discovered that liquid 4 He below the A-point can 
flow almost frictionless. Rapitza named this behavior 
superfluidity [TJ [2] . Many of the properties of superfiuid 
helium also appear in the gaseous Bose-Einstein conden- 
sates (BECs). In contrast to liquid helium, where the 
interatomic interaction is too strong to investigate the 
microscopic properties of superfluidity, the interactions 
in gaseous BECs are much weaker. The study of su- 
perfiuid flow in dilute BECs can therefore deepen our 
understanding of superfluidity. 

The observation of vortices in a gaseous BEC [3] H] 
and the demonstration of the persistent flow of a BEC 
in a toroidal trap [5] give a striking demonstration of 
superfluidity. Evidence of the breakdown of superfluidity 
in a BEC is obtained in an experiment by Raman et al. , 
in which a blue detuned laser beam is moved through the 
BEC at different velocities [6]. Strong heating is observed 
only above a critical velocity v ci which is found to be 
v c ~ 0.25c, with c the local speed of sound at the peak 
density of the BEC. 

In this Letter the flow of a BEC through a thermal 
cloud in a harmonic potential is studied by exciting a 
dipole oscillation of the BEC, whereas the thermal cloud 
initially remains stationary. In the hydrodynamic regime 
this out-of-phase mode of the trapped Bose gas is the 
analog of the usual second sound mode in bulk superfiuid 
helium [7J. For this second sound dipole mode we study 
for the first time its frequency and damping rate from the 
collisionless to the hydrodynamic regime. In contrast to 
liquid helium, our analysis allows for a direct measure- 
ment of the position of the superfiuid component (con- 
densed atoms) with respect to the normal fluid (thermal 
atoms), which allows for an unequivocal determination 
of the second sound dipole mode. A dipole oscillation of 
the center-of-mass of the whole system is undamped in a 
harmonic potential and this in-phase dipole oscillation is 
often used to measure the trap frequency. 



The thermal cloud can also be tuned into the collision- 
less regime, in which the mean free path of the thermal 
atoms is larger than the axial size of the cloud. As we will 
show, the damping of collective excitations in a collision- 
less partly condensed BEC is primarily caused by Landau 
damping, i.e. mean-field interactions mediate the transfer 
of energy from the condensate to the thermal cloud, lead- 
ing to the damping of collective modes. Landau damping 
was first discussed by Landau in the context of the damp- 
ing of plasma oscillation and plays a key role in a broad 
variety of fields, for instance the damping of phonons in 
metals, the damping of quarks and gluons in quark-gluon 
plasmas and the anomalous skin effect in metals. 

Previously, some experiments have been performed in 
which the BEC and the thermal cloud move with respect 
to each other and either the BEC [5] or the thermal cloud 
|9J is initially at rest. In the experiment presented in 
Ref. [9] the thermal cloud is in the collisionless regime. In 
a pioneering experiment presented in Ref. [8 the thermal 
cloud is in the crossover region toward the hydrodynamic 
regime, but only one measurement series is presented. It 
therefore does not provide a study of the dependence of 
the superfiuid flow on for instance the temperature or 
the collision rate. In contrast, the number of theoretical 
studies of this subject greatly outnumbers the available 
experimental work [7J [TO] EEU H EES] • 

The experimental setup, as described in detail in 
Ref. [TJ], is capable of creating Bose-Einstein conden- 
sates containing up to 3 • 10 8 sodium atoms, limited by 
three-body decay. In this setup we have reached the axi- 
ally hydrodynamic regime in the thermal cloud above the 
transition temperature T c . The axially hydrodynamic 
regime is reached by evaporatively cooling atoms in an 
axially strongly decompressed trap with an aspect ratio 
of up to 1 : 65 [IB] • The radial trap frequency is given 
by uj ra d/(2n) — 95.56 Hz. Here, we cool atoms to tem- 
peratures below T c for three different values of the axial 
confinement. 
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We introduce a measure for the hydrodynamicity of 
the thermal cloud in the axial direction 7 = 722 /wax, 
where the collision rate 722 = n e s a w roi is the average 
number of collisions in the thermal cloud. Here, v re \ 
is the relative velocity, a = 87ra 2 is the isotropic cross- 
section of two bosons with s-wave scattering length a and 
?^cff = / drn 2 x (r)/ J dfn ex (f) is the effective density of 
the trapped thermal atoms with density n ox . For our pa- 
rameters, 722 — 90s -1 for the highest number of atoms 
and weakest axial confinement, which corresponds to a 
hydrodynamicity of 7 < 10. Note that the thermal cloud 
becomes less hydrodynamic for stronger axial confine- 
ment due to two effects: the increased three-body losses 
and the radial expansion. We have observed the crossover 
from the collisionlcss regime to the hydrodynamic regime 
by studying the heat conduction in a thermal cloud above 
T c [IB]. Furthermore, it is noted in Ref. [TT] that the ther- 
mal cloud is deeper in the hydrodynamic regime when a 
BEC forms due to collisions with condensed atoms. As 
a result the thermal cloud is even more hydrodynamic 
below T c than it is above the transition temperature. 




FIG. 1: Schematic representation of the excitation procedure 
of the second sound dipole oscillation. See text for more de- 
tails. 

The procedure used for the excitation of a second 
sound dipole mode, schematically shown in Fig.[l] is con- 
ducted as follows: initially, the center of the magnetic 
trap (MT) at z — z (Fig. [Ik) is displaced adiabatically 
in the axial direction over slightly more than the length 
of the BEC (Fig.[T|D). Next, a focused blue-detuned laser 
beam is aligned perpendicular to the axial axis of the sys- 
tem at position with z$ < < z±, where it does not 
overlap with the BEC (Fig. [II). This beam is detuned 
25 nm below the 23 Na D2-transition and its intensity 
is gradually ramped up to roughly 10 4 mW/cm 2 in 100 
ms. This corresponds to a maximum potential of roughly 
Vdip — 1.1 x/i, where \i is the chemical potential of the ini- 
tial cloud. The intensity is continuously measured using 
a photo diode in order to stabilize the intensity using a 
feed-back circuit controlling the efficiency of an acousto- 



optical modulator, which deflects the light used for the 
dipole beam. Next, the center of the trap is displaced adi- 
abatically back to its original position z = zq (Fig. [TJi). 
Since /j, < Vdi p < k^T, with &b the Boltzmann constant 
and T the temperature, only the thermal cloud can pass 
the dipole potential and returns to z = zq, while the con- 
densate remains at z = z\. Finally, the dipole beam is 
turned off and the BEC will undergo a dipole oscillation, 
while the thermal cloud is initially at rest (Figs. [XJs and 
[Tj 7 ). The distance Az = z& — zq determines the maximum 
velocity v = w ax Az of the center-of-mass of the conden- 
sate. The final position of the center of the MT is used 
to tune Az and thus v. 

We have excited the cloud for various v in order to ob- 
serve, if there is a threshold velocity v c above which the 
out-of-phase oscillation is strongly damped. The damp- 
ing rate for different values of v is shown in Fig. [2j For 
v > c we observe an increase of the damping rate of 
the out-of-phase oscillation for increasing v. The strong 
damping in this regime is accompanied by a large reduc- 
tion of the condensate fraction, which eventually leads to 
the complete depletion of the BEC. For smaller v (but 
still above v c ) we observe two distinct damping rates as 
long as the relative velocity is sufficiently reduced below 
v c before the condensate is fully depleted: first we observe 
a high damping rate, but as soon as the relative veloc- 
ity is below v c a much lower damping rate is observed. 
For excitations with v < v c we observe a low damping 
rate, independent of v within the experimental accuracy. 
The threshold value depends on our trap geometry and 
is determined at v c /c « 0.59 ± 0.05 for the weakest axial 
confinement. For the strongest confinement used we find 
v c /c w 0.8 ± 0.1. The observation of a threshold veloc- 
ity indicates the onset of a strong damping mechanism, 
which we consider as the breakdown of superfluidity. The 
experiments described in the remainder of this Letter are 
conducted in the regime of low damping (v < v c ). 
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FIG. 2: The damping rate of the out-of-phase oscillation as 
a function of the initial velocity v/c, where c is the central 
speed of sound determined using the measured BEC density. 
The uncertainty in the measured damping rate and in c is 
larger for high v due to the depletion of the condensate. 

A measurement series contains roughly 100 shots at 
various hold time r after the dipole beam is suddenly 
turned off. At the end of the hold time an absorption 
image is taken after the cloud is allowed to expand dur- 
ing a 85 ms time-of- flight. During a measurement series 
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the situation before the dipole beam is turned off is mon- 
itored a few times to check for the absence of condensed 
atoms in the center of the trap. The thermal cloud is eas- 
ily distinguished from the BEC in the absorption images 
due to the distinct density profiles. The final images are 
fitted to a bimodal distribution. This distribution is the 
sum of a Maxwell-Bose distribution modeling the ther- 
mal cloud and a Thomas-Fermi (TF) distribution model- 
ing the BEC. The fit yields ten parameters; the axial and 
radial positions, the axial and radial sizes and the optical 
densities of both components. The axial position of both 
components is used to observe the dipole oscillations, the 
total absorption is used to determine the column density 
of both components, and the radial size of the thermal 
cloud is used to determine the temperature of the sys- 
tem. Furthermore, the aspect ratio of the condensate is 
used to detect the quadrupole oscillation induced by the 
excitation procedure, which turns out to be small. This 
mode does not couple to the dipole mode and is therefore 
not of interest to us here. 
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FIG. 3: The in-phase and out-of-phase dipole oscillation of a 
system characterized by T — 0.39/xK and fi/h = 2.6 kHz for 
Wax/ (27r) = 7.78 Hz. In Figs, (a) and (b) the axial position 
z cx of the thermal cloud and zbbc of the BEC is plotted as 
a function of the hold time r. In Fig. (c) the motion of the 
out-of-phase dipole oscillation is isolated by plotting Az as 
a function of r. The solid lines are the result of a fit to the 
data. Due to the destructive imaging scheme used each point 
represents the position of a newly prepared cloud. 



A series is measured for three different values of the 
axial confinement and for two different temperatures. To 
compare the results under these different conditions, the 
temperature is adjusted in such a way that for a given 
temperature the initial chemical potential is roughly the 
same in all confinements. The axial position of the 
thermal cloud z cx and the BEC zbec are plotted as 
a function of r in Fig. [3] Furthermore, the difference 
Az = zbec — z cx is plotted to isolate the out-of-phase 
dipole oscillation from the in-phase oscillation. The ther- 
mal cloud, initially at rest, starts to oscillate when the 



BEC moves through it and eventually both the BEC and 
the thermal cloud oscillate simultaneously and in-phase, 
as can be seen in Fig. [3j We refer to this final motion 
as the in-phase dipole oscillation. We have checked that 
the thermal cloud remains at rest as long as the dipole 
beam is present. 

The positions z ex and zbec are simultaneously fitted 
to a combination of two exponentially damped sinusoidal 
functions. The fit of the data yields the damping rate 
f"ip (Toop) and the frequency u>i p (w op) of the in-phase 
(out-of-phase) dipole oscillation. The damping rate of 
the in-phase dipole oscillation does not deviate from zero 
within the experimental uncertainty, as is expected for 
a harmonic trap. The axial trap frequency determined 
from the in-phase oscillation frequency gives the same 
result as the trap frequency determination from a mu- 
tual center-of-mass oscillation (tJi p = w ax ). The quanti- 
ties r oop and w op are plotted as a function of the trap 
frequency in Fig. [4] 
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FIG. 4: The damping rate r oop (a), the normalized damp- 
ing rate r oop /w ax (b), and the normalized oscillation fre- 
quency Lj oop /ijj ax (c) as a function of the axial trap frequency 
w ax /(27r). The diamonds (triangles) indicates the high (low) 
temperature condition. 

In Fig. [4]i one can see that r oop decreases for decreas- 
ing axial trap frequency. However, r oop /aj ax increases for 
decreasing axial trap frequency as one can see in Fig.[3]3. 
This indicates that an extra damping mechanism, which 
does not depend on the trap frequency, plays a signifi- 
cant role as well, especially for the highest T and thus 
the most hydrodynamic clouds. 

The normalized second sound oscillation frequency 
^oop/w ax is shown in Fig. For all conditions uj oop < 
w ax , where under the most hydrodynamic conditions 



0.94w ax . Furthermore, the largest deviations of 
/w ax from 1 coincides with the largest r oop and is 
found for the most hydrodynamic clouds. The shift of 
with respect to w ax , up to 6% in our experiment, 



^oop 
^oop 



^oop 



is larger than the shift of 4.5% observed in Ref. [5] and 
similar to the 6% shift reported in Ref. |9J, although the 
reduced temperature T/T c is lower in our experiments. 
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The mean-field potential of the thermal cloud causes 
the effective axial trapping potential to be slightly 
smaller. This effect reduces the effective trap frequency 
with roughly 2% for the highest thermal density in these 
experiments and cannot account for the shifts of up to 
6%. As a consequence, the measured frequency shift and 
damping rate of the out-of-phase dipole mode both re- 
flect that collisional effects play a role. Since uj oop and 
u> ax are determined simultaneously and with a high ac- 
curacy the frequency shift gives an accurate measure for 
these effects. 

The measured damping rates are compared in Fig. [5] 
to Landau damping. Landau damping in a homoge- 
neous, weakly interacting Bose gas has been described 
in the collisionless limit by several authors [TTJ [TS] . Lan- 
dau damping increases with temperature, because of the 
larger number of particle-hole pairs available at thermal 
equilibrium |19J . In our case, due to the confinement 
both the BEC and the thermal cloud have an inhomo- 
geneous density distribution. Landau damping in such 
a system is determined by the condensate boundary re- 
gion, and the result for the damping rate is different 
from that in a spatially homogeneous gas (2Qj EH 122] , 
The damping rate for an inhomogeneous gas is given 
by Ll = CLljj = C Sirk^T a u> / (8Hc) , where the speed 
of sound is given by c = y/ fi/m, with m the atomic 
mass [22]. Furthermore, u> is the frequency of the excita- 
tion, is the Landau damping rate for a homogeneous 
gas and C = 12cto/(7t 3 / 2 fee Vmhujlog [2c 2 m/(hu>)] ) is 
a numerical coefficient to incorporate the inhomogcneity. 
Note that in our experiment r L is roughly independent 
of the speed of sound, since fj, and therefore the central 
density is roughly constant in the different axial confine- 
ments. 

Fig. [5] shows that we have measured the damping rate 
of the second sound dipole mode in a range from close to 
the collisionless regime up to the hydrodynamic regime 
with a hydrodynamicity close to 10. The ratio r oop /rL 
increases strongly as a function of the hydrodynamicity 
7 and this may serve as a stringent test for theoretical 
models, which attempt to describe the damping of hydro- 
dynamic excitations in the Landau two-fluid (normal and 
superfluid) model. It shows that collisional processes, 
which have not been incorporated in the descriptions of 
Landau damping discussed above, play an important role 
in damping the second sound dipole mode. 

In conclusion, we successfully excited a second sound 
dipole mode in a trapped Bose-condensed gas at finite 
temperature. The damping rate and the oscillation fre- 
quency of the second sound dipole mode is measured as 
a function of the axial confinement and the temperature 
with a high accuracy. We find that Landau damping, 
when the inhomogeneous density of the cloud is taken 
into account, shows good agreement to the measured 
damping rates close to the collisionless regime. In the 
hydrodynamic regime the measured damping rates are 
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FIG. 5: The ratio of the measured damping rate and Landau 
damping r oop /rL as a function of the hydrodynamicity pa- 
rameter 7. The data point represented by the triangle is the 
result of a measurement taken from Ref. [S] . 



larger than the rates based on Landau damping. We hope 
that the measurements presented here encourages calcu- 
lations based on models that take collionsional damping 
into account and can provide a fruitful testing ground for 
such models. 
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